On the Origin of Divergence Errors in Mhd Simulations and Consequences for Numerical Schemes
نویسنده
چکیده
A. This paper investigates the origin of divergence errors in MHD simulations. For that purpose, we introduce the concept of discrete involutions for discretized conservation laws. This is done in analogue to the concept of involutions for hyperbolic conservation laws, introduced by Dafermos. By exploring the connection between discrete involutions and resonance, especially for constrained transport like MHD, we identify the lack of positive central viscosity and the assumption of one-dimensional physics in the calculation of intercell fluxes as the main sources of divergence errors. As an example for the consequences for numerical schemes, we give a hint how to modify Roe-type schemes in order to decrease the divergence errors considerably and, thus, stabilize the scheme.
منابع مشابه
Increasing the accuracy in locally divergence-preserving finite volume schemes for MHD
It is of utmost interest to control the divergence of the magnetic flux in simulations of the ideal magnetohydrodynamic equations since, in general, divergence errors tend to accumulate and render the schemes unstable. This paper presents a higher-order extension of the locally divergence-preserving procedure developed in Torrilhon [M. Torrilhon, Locally divergence-preserving upwind finite volu...
متن کاملStability Analysis and Error Estimates of an Exactly Divergence-Free Method for the Magnetic Induction Equations∗
In this paper, we consider an exactly divergence-free scheme to solve the magnetic induction equations. This problem is motivated by the numerical simulations of ideal magnetohydrodynamic (MHD) equations, a nonlinear hyperbolic system with a divergence-free condition on the magnetic field. Computational methods without satisfying such condition may lead to numerical instability. One class of me...
متن کاملA constrained-gradient method to control divergence errors in numerical MHD
In numerical magnetohydrodynamics (MHD), a major challenge is maintaining ∇ · B = 0. Constrained transport (CT) schemes achieve this but have been restricted to specific methods. For more general (meshless, moving-mesh, ALE) methods, ‘divergence-cleaning’ schemes reduce the ∇ · B errors; however they can still be significant and can lead to systematic errors which converge away slowly. We propo...
متن کاملLocally Divergence-preserving Upwind Finite Volume Schemes for Magnetohydrodynamic Equations
A main issue in nonstationary, compressible magnetohydrodynamic (MHD) simulations is controlling the divergence of the magnetic flux. This paper presents a general procedure showing how to modify the intercell fluxes in a conservative MHD finite volume code such that the scheme becomes locally divergence preserving. That is, a certain discrete divergence operator vanishes exactly during the ent...
متن کاملAIAA 2001–2623 Multi-Dimensional Upwind Constrained Transport on Unstructured Grids for ‘Shallow Water’ Magnetohydrodynamics
Novel Multi-dimensional Upwind Constrained Transport (MUCT) schemes on un-structured triangular grids are described. Constrained Transport (CT) discretizations conserve the divergence-free nature of divergence-free vector fields on the discrete level. Multi-dimensional Upwind (MU) schemes generalize the concept of dimensionally split upwind schemes for hyperbolic systems to truly multidimension...
متن کامل